Frozen local hole approximation 
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The frozen local hole approximation (FLHA) is an adiabatic approximation which is aimed to 
simplify the correlation calculations of valence and conduction bands of solids and polymers or, 
more generally, of the ionization potentials and electron affinities of any large systems. Within 
this approximation correlated local hole states (CLHSs) are explicitely generated by correlating 
local Hartree-Fock (HF) hole states, i. e. (N-l )-particle determinants in which the electron has 
been removed from a local occupied orbital. The hole orbital and its occupancy is kept frozen 
during these correlation calculations, implying a rather stringent configuration selection. Effective 
Hamilton matrix elements are then evaluated with the above CLHSs; diagonalization finally yields 
the desired correlation corrections for the cationic hole states. We compare and analyze the results of 
the FLHA with the results of a full MRCI(SD) (multi-reference configuration interaction with single 
and double excitations) calculation for two prototype model systems, (H.2) n ladders and H-(Be) n -H 
chains. Excellent numerical agreement between the two approaches is found. Comparing the FLHA 
with a full correlation treatment in the framework of quasi-degenerate variational perturbation 
theory reveals that the leading contributions in the two approaches are identical. In the same 
way it could be shown that replacing a correlation calculation around a frozen local hole is in 
fact equivalent, up to first order, to perform a much less demanding SCF (self-consistent field) 
calculation around the hole to relax all other occupied orbitals. Thus, both, the FLHA and the 
above SCF approximation, are well-justified and provide a very promising and efficient alternative 
to fully correlated wavefunction-based treatments of the valence and conduction bands in extended 
systems. 

PACS numbers: 71.10.-w,71.15.-m,71.20.-b 



I. INTRODUCTION 

In the last decade, increasing interest in wavefunction- 
based correlation methods for excited states in extended 
systems can be observed HSHSlSHSSEil- 
This is because one can rely on the very sophisticated nu- 
merical methods for molecules which are well-established 
in quantum chemistry. These methods arc conceptually 
very clear and allow for systematic improvement: based 
on self-consistent field calculations, electron correlation 
can be included successively by considering determinants 
of increasing excitation order in configuration interaction 
(CI) or related procedures. They yield approximate cor- 
related wavefunctions from which all quantities of inter- 
est can be derived. 

Of course, wavefunction-based methods are computa- 
tionally very demanding. In order to arrive at manage- 
able schemes for the description of solids or other ex- 
tended systems, local orbitals have to be introduced and 
the local character of electron correlation must rigor- 
ously be exploited. Among these local approaches \V\\ 
one can distinguish local MP2 (M0ller-Plesset perturba- 
tion theory up to 2 nd order) and coupled electron pai r 
or coupled cluster schemes [H OS Q [H 03, El, 
methods which are based on a Green's function formal- 
ism pi, ^3 an d techniques which use local Hamiltonian 
matrix elements El La La III or other partitioning concepts 
H3, 123, 123, ISjnSa |SE [53, |2l, |2^|. Whereas the 
first-mentioned methods work in infinite systems and aim 
to truncate the excitation space by an appropriate con- 
figuration selection scheme, the effective Hamiltonian ap- 



proach as well as the other partitioning schemes assemble 
the correlation effects from finite subsystems (embedded 
fragments) which are finally transfcrcd to the infinite sys- 
tem. Most of these methods focus on ground state prop- 
erties, only the Green's function and local Hamiltonian 
approaches are explicitely designed for excited states. In 
particular, the local Hamiltonian approach has been used 
in the past quite successfully to describe the correlation 
effects on the band structure of covalent solids 0, 0, 
and polymers . Yet, this approach is quite involved and 
a more simpler, approximative way of computing the lo- 
cal Hamiltonian has been proposed recently by our lab 

M- 

In this work we want to analyze the guiding approxima- 
tions introduced in that recent approach, and therefore 
address the question to which extent it is possible to sim- 
plify the correlation calculations by focusing on "frozen" 
local hole configurations although, in reality, the electron 
hole is usually delocalized over the entire system. Intu- 
itively, this frozen local hole approximation, is guided by 
the idea that the shape of the correlation hole which is 
carried along by a traveling electron is essentially invari- 
ant and is thus following the electron hole adiabatically. 

A related approach is the frozen core hole approxi- 
mation [3lJ for the calculation of core-hole states where 
the core hole orbital is frozen and all other orbitals 
are allowed to relax during an Hartree-Fock calculation. 
Such a local view on core holes in bulk materials is also 
adopted in other wavefunction-based investigations on 
excited core hole states [32, While it seems quite 
natural to describe core holes in that way, although in 
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solids even core orbitals are Bloch waves, it is much less 
obvious that a frozen local hole treatment should also be 
possible for delocalized valence hole states. 

In the following, we will outline the theoretical back- 
ground of the frozen local hole approximation (for valence 
holes) and then present and analyze the numerical results 
obtained for two model systems, namely (H2)„ ladders 
and H-(Be) n -H chains. These two examples have been 
chosen to represent the two cases of prevailing van der 
Waals and predominantly covalent binding. Finally, an 
analytical, perturbative analysis of the frozen local hole 
approximation is given. 

II. THEORETICAL BACKGROUND 

Starting point for our correlation calculations is the 
Hartree-Fock (HF) ground state Slater determinant \&q) 
and energy Eq of the iV-particle system. The canonical 
HF orbitals \v) with orbital energies e„ are divided into 
the N/2 energetically lowest-lying orbitals which are oc- 
cupied in |$Q } and the unoccupied, virtual orbitals. Fur- 
thermore, the occupied orbitals are subdivided into core 
and valence orbitals depending on whether they are kept 
frozen or are active during the subsequent correlation 
calculation. Within these orbital groups we will switch 
between delocalized canonical orbitals and localized or- 
bitals, {|^)} and {|a)}, using unitary transformations, 
such as 

\a)=J2\»)U»a (1) 

v 

for the occupied orbitals. In order to find the correlated 
cationic (A r -l)-particle states of the system in mind, one 
first constructs the so-called reference states l^f^ 1 ) of 
the system by removing one electron out of one of the 
valence orbitals = Cj|<I>^). The so generated 

canonical hole states l^^" 1 ) and local hole states 
are connected via 

K- 1 } = E l^ 1 )^ 1 ):, = E l^a-'Wua (2) 
a a 

Since we are only dealing with cationic states here we 
will omit the upper index (N-l) in the following. The 
canonical hole states are the eigenstates of the CI matrix 
in the reference space. They describe the hole states of 
the system according to Koopmans' theorem |34| . 

Including correlation effects one arrives at correlated 
wave functions with energies E v . In analogy to 

Eq. J5J we then define so-called correlated local hole states 
|¥„) (see Fig. [TJ 

l*.):=El*")^ ^ 
and an effective Hamilton matrix if cff = (H^b) w i^h 

ff*:=(fa|^l$i) = ^Wi ■ (4) 



(I**)} ^ HF states 

i 

MRCI j FLHA 

T 
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FIG. 1: Relation between the delocalized hole state on the 
Hartree-Fock level [<&„) and the corresponding correlated 
states \9 V ), and their localized counterparts, |$ a ) and \^ a ), 
respectively. Within the frozen local hole approximation 
(FLHA) approximate correlated local hole states |\I/ a ) are 
constructed without explicit reference to the true correlated 
states 



The effective Hamiltonian H is constructed such that 
diagonalizing its local matrix representation H_ cS pre- 
cisely yields the desired correlated energies E v of the 
cationic states. Thus, provided one knows the correlated 
local hole states |* Q ) defined in Eq. (J3J, one can find 
the cationic energies by simply solving a small eigenvalue 
problem. 

Within our frozen local hole approach we do not follow 
the three-step process depicted in Fig.^but directly gen- 
erate approximate correlated local hole states (CLHSs) 
\Sfr a ) by performing separate correlation calculations for 
each reference state |<E> a ) during which the hole in the 
localized orbital \a) is kept frozen. This implies a con- 
figuration selection where only those configurations are 
used in which the hole in \a) is maintained. In the case 
of an CI(SD) (CI with single and double excitations) cal- 
culation the approximate localized hole states \^f a ) take 
on the following form: 

|*a) = «a|$a)+E<,J^,x)+E<i'l $ «£x') " ( 5 ) 

x,v x ,x' 
v,v' 

The indices x,x' € {a,b,b, ...} run over all the remain- 
ing valence electrons, with a being the electron with op- 
posite spin to the removed electron a. v, v' denote the 
electrons in the virtual orbitals. Thus, e. g., a ) is 
the 2hlp configuration with no electrons remaining in 
the spatial orbital associated with a whereas |<I>^ b ) con- 
tains holes in two different spatial orbitals. Obviously, 
the configuration spaces for the different hole states \^> a ) 
are overlapping. Excited configurations with holes in two 
or more valence orbitals |o), \b), ... are, in fact, present in 
the correlation calculations of each reference state |$ Q ), 
. . . (see Fig. El . As we will see in detail later, this 
fact does not harm the procedure. On the contrary, it is 
important to build up the correlation effects in the final, 
delocalized hole states properly through a suitable mix- 
ing of the configurations in the final diagonalization step 
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FIG. 2: Geometry of the (H2) n ladders used in this paper. 
The H2 bond length is fixed at its molecular value of R = 
0.7417 A; the distance between the H2 units is varied during 
the calculations. 

of the effective Hamiltonian Q , 

|#*)=X>a(>)l f a) where (6) 

a 

J2Hlb^n=E u ^S ab X b (u) . (7) 

b b 

Note, that the individual approximated CLHSs \*& a ) are 
not orthonormal with respect to each other. Thus, one 
has to solve a generalized eigenvalue problem here (with 
S a b = (^a\^b) being the overlap matrix) in order to ar- 
rive at the approximate canonical hole states and 
the approximate hole state energies E v . 

III. APPLICATIONS 

As first applications of the frozen local hole approxima- 
tion (FLHA) we choose two simple model systems, (H2)„ 
ladders (see Fig. [5J and linear H-(Be)„-H chains, for 
which we compare and analyze the results of the approx- 
imation with the results of the corresponding complete 
multi-reference CI(SD) (MRCI(SD)) calculation. All re- 
sults are obtained by using the MOLPRO program pack- 
age in particular the MRCI(SD) option IH l3c| . 
For the calculations sp cc-pVDZ basis sets by Dunning 
[3^| are used for H and Be except for the terminating 
H atoms in the Be chains which are described by a re- 
duced s cc-pVDZ basis. The Be-H distance has been 
set to 1.37 A which is the equilibrium distance found in 
all neutral chains of length rt=4-10. The computation 
of the Hartree-Fock ground state of the neutral system 
\&q) yields the canonical Hartree-Fock orbitals. The set 
of valence orbitals is subse quen tly localized by means of 
the Foster-Boys procedure [4(| and is used for the con- 
struction of the localized hole states {|$ a )}. Each |$ a ) is 
then separately correlated in a CI(SD) calculation with 
the above-described configuration selection, where only 
those single and double excitations are allowed with a 
hole in the localized orbital \a). In order to achieve this 
configuration restriction within the MOLPRO program, 
one has to declare all valence and unoccupied orbitals as 
active orbitals. The maximum number of active orbitals 
is fixed to 32 in the MOLPRO code, which only allowed 
us to handle relatively small chains and basis sets chosen 
here. Nevertheless, the effects of the approximation can 
be studied very well in these model systems. 

The results for the (H2)^~ and (H2);j~ ladders are shown 
in Figs. and 0] Here, the total energies of the hole 
states on the HF, the MRCI(SD) and the FLHA level 
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FIG. 3: Total Hartree-Fock, MRCI(SD) and FLHA energies 
of (H2)J as a function of the H2-H2 distance. 
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FIG. 4: Total Hartree-Fock, MRCI(SD) and FLHA energies 
of (H2)^" as a function of the H2-H2 distance. 



are depicted in dependence of the distance between the 
H2 units. In accordance with the number of valence or- 
bitals, and thus possible hole states, we find three and 
five low- lying cationic states of (H 2 )J and (H 2 )^, respec- 
tively. In their lowest states the cationic (H2)„ ladders 
are stable. Compared to the HF data, very pronounced 
correlation effects are observed: In the case of the (H2) 3 
chain the lowest-lying state is lowered by about 2.9 eV 
at its equilibrium distance, which is slightly shifted from 
i? oq = 1.7 to 1.6 A, and the potential well is deepened 
from 1.6 to 1.9 eV. Also for the high-lying states of (H2)^ 
and the states of (H 2 )^" states we find an overall energeti- 
cal lowering of about 3 eV through correlation. The most 
important finding in the present context, however, is that 
all the described effects are very well accounted for by the 
frozen local hole approximation. The FLHA curves fol- 
low the MRCI reference data very closely with the only 
noticeable slight deviations at the outer left edge of the 
potential curves. 

In order to make sure that this excellent agreement 
is not only an artefact of the relatively weak van der 
Waals binding between the H2 units, we choose the 
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FIG. 5: Total Hartree-Fock, MRCI(SD) and FLHA energies 
for Be3H^ in dependence of the Be-Be distance. 



predominantly covalently bound H— (Be) n — H chains as 
second application. For all systems considered, H 
(Be) 3 -H through H-(Be) 5 -H, the FLHA and the "exact" 
MRCI(SD) data agree very well. In the following, we only 
want to discuss in detail the results for the smallest sys- 
tem Be3rlJ shown in Fig. The calculated correlated 
potential energy curves for the two lowest cationic states 
exhibit shallow minima at Be-Be distances of R cq = 2.3 
and 2.4 A. For the most stable cationic state of Be3H2 
the correlation energies vary from about 1.1 eV at the 
minima to about 1.5 eV for small and large Be-Be dis- 
tances. 

The dominating configurations of the two lowest 
cationic states are the two reference configurations |<& a ) 
and |4>b) which result from taking one electron out of 
one of the highest occupied molecular orbitals (HOMOs) 
of the HF determinant. The corresponding localized or- 
bitals are centered at the Be-Be bonds (see top row of 
Fig.EJ. Within the FLHA two configuration spaces I and 
II are constructed, built upon the two reference config- 
urations |$ Q ) and l^b). They are shown schematically 
in Fig. |SJ Clearly, the two configuration spaces strongly 
overlap. For this simple example with two reference con- 
figurations only, all mixed (single-excited) 2hlp configu- 
rations |$^), |$a fc ) and and all (doubly-excited) 

3h2p configurations \$ v a ,v xx i) and \$l' x x ,} fall into the 
overlap region. Nevertheless, very important configura- 
tions, especially the reference configurations themselves, 
lie in only one of the spaces. That is why this very sim- 
ple example is already suited to show the main features 
of the approximation. Strictly speaking, the single de- 
terminants \& v z) only belong to space I (as the b ) 
only belong to space II). But since MOLPRO uses spin- 
adapted excitations these two types of single determi- 
nants always come in singlet and triplet- like linear combi- 
nations (|$^ 5 ) ± and are thus both attributed 
to the overlap region I fl II. 

In the present example, no substantial saving in the 
computational cost can be gained from using the FLHA, 




FIG. 6: Sketch of the overlapping configuration spaces I and 
II for a system with two valence orbitals and two hole states. 
The nomenclature is the same as in Eq. @. 



since we substitute the single MRCI(SD) calculation by 
two CI calculations with configuration spaces of almost 
the same size as the original one. But this will change 
dramatically with bigger systems because a given config- 
uration can be found in at most three different config- 
uration spaces. E. g., the double excited configuration 
Italic) would only be part of the configuration spaces 
generated from the reference configurations |$ a ), |$&) 
and I'&c)- Once suitable cut-off criteria are introduced, 
the configuration spaces of each reference configuration 
remain finite no matter how large the systems becomes; 
and since the number of these spaces only grow linearly 
with system size the FLHA directly leads to an O(N) 
method. 

Coming back to the discussion of the Be3H2 results, 
only a few configurations with noticeable weights (> 
0.05) are found in the approximated correlated local hole 
states |* Q ) and (Eq. ©). The same holds true 

for the "full" canonical MRCI wave functions, l*^ 1101 ) 
and |^2 RCI ). In order to enable an analysis of the indi- 
vidual correlation contributions, we also localize the vir- 
tual orbitals, although this is by no means a prerequisite 
of our FLHA. All relevant configurations contain excita- 
tions into one of the four pair-wise degenerated virtual 
orbitals shown in Fig. [7| The first pair, named c* and 
c*, is localized at the central Be atom, the other pair, 
a* and b* , represents anti-bonding orbitals on the two 
Be-Be bonds. 

In Table [I] we summarize the CI coefficients ai and 
[3i of the most relevant contributions to the approxi- 
mate CLHSs y$> a ) an d |^h)j as well as the CI coefficients 
(aj ± j3i)/y/2 of the approximated wave functions , F„ of 
the hole states 1 and 2 as resulting from the diagonal- 
ization of the effective Hamiltonian, a (2 x 2) matrix, 
here. The latter CI coefficients are compared with those 
of the "exact" MRCI wave functions, aff cl and aff 01 . 
For all spin-adapted single determinants considered here: 
the reference configurations, the single and double exci- 
tations into the anti-bonding orbitals a* and b* and the 
3h2p configurations involving the orbitals c* and c*, an 
excellent agreement is found between the approximate 
and the "exact" CI coefficients of the holes states with 
deviations well below 10% in most of the cases. In con- 
trast to our initial expectation, no cancellation of artifi- 
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e=-0.349 





c 8=0.244 c 




a* e=0.423 b 

FIG. 7: Localized highest occupied orbitals (HOMOs) of 
BesHj together with those localized virtual orbitals which 
are involved in the most relevant configurations of the ap- 
proximate correlated local hole states \^ a ) an d |^f>)- The 
Hartree-Fock orbital energies e of each energetically degener- 
ate pair of orbitals are given as well. 



where Q = J^i \&i) (®i I is the projector onto the or- 
thogonal complement of the model space and the con- 
figurations run over all single, double and higher 
excitations one can form from the local model space con- 
figurations {|$ Q )}. 

In first order perturbation theory the wave operator Q 
is given by 1^3 



n - X>c>(*c| 

c 

+ E E i*i 



(10) 



z\\ab) 



c a<b,v 



£b 



We take the zeroth-order Hamiltonian H° to consist of 
the diagonal terms of the Fock operator only, i. e., Hy = 
SijEi with Si = Fa. The perturbation then contains both, 
the off-diagonal terms F„ = (<pi\F\<pj) , i ^ j of the 
Fock operator F and the usual two-electron contributions 
(ij\\kl) — X)a CC (* a llj a )- Here and in the following, indices 
a, b, ... are understood to run over occupied localized 
orbitals while indices v, w, ... run over virtual orbitals. 
For arbitrary orbitals indices i, j, . . . are used. The 12,12 
notation 



cial, orbital-relaxation-based contributions in the FLHA 
wave functions occur. All configurations which are im- 
portant in the CLHSs also contribute significantly to (at 
least) one of the final hole states |^i) and | \I/ 2) ■ Appar- 
ently, the ultimate diagonalization step just mixes the 
individual contributions in such a way that the delocal- 
ized nature of the hole states is recovered, precisely as 
one would expect for an adiabatic approximation. 



(ij\\M) = (mid) - mm with (11) 

m\kl) =U P?(1M(2) r£ p*(l) W (2) dld2 

is employed throughout. The |<^,-) are the (localized) 
Hartree-Fock orbitals of the neutral A-electron system. 
Using that form of the wave operator as variational 
ansatz for the eigenstates of the hole system, i. e. 



IV. PERTURBATIVE ANALYSIS |*„) = Q[ ^ A (i/)|$ ) ] = E X "( v ) > ( 12 ) 

a a 

To understand, why the frozen local hole approxima- 
tion (FLHA) works so well, we switch from the MRCI leads to the following effective secular matrix 
level of theory to quasi- deg enerate variational perturba- 
tion theory (QDVPT) 41] and try to find the leading H°f := (fi[$ e ] \H\ &[® f ] ) (13) 
contributions in both, the "exact" hole state and the ap- 
proximate hole states according to the FLHA. Within = — Ffe (14) 
QDVPT the correlated wave functions \^ v ) of a system , H r v / f\\ u\ 

... ' / . VV e\ab,f \ v J\\ ao / 

are written as '—r 1 

l*"> = A l*?> ^ + "e We*vW\ve) 

where = P\^u) are the projections of the corre- a<b,v 

lated wave functions onto the space M spanned by _|_ \^ W e v ab \^ d j (ab\\ve) (wf\\cd) 

the local Hartree-Fock hole configurations {|$ a )}, and Cl a <b,vc<d,w 

is an operator which acts on the "model space" M. and with 

provides the full, correlated wave functions \^ v )- It is ir/isti \ 

assumed that the \^ v ) are those eigenstates of the hole W e] v ab j = - — V e| 1 ahl — - (15) 

system which are dominated by the single determinants ^ £v £a £b £ 5> 

forming the model space and, in particular, that the pro- (<I>^ b |i?|$/) v 

jections \^{?) remain linear independent. In fact, e,ab\f • i £ ^ _ £a _ £fe _|_ £ ^ ~ <• f\ab,e) 

n = P + QCl with QCl = Y |$/)fi/ a ($a| (9) w v w = (Kb\H\®cd 



e,ab\cdj - _ ^ _ £b + £ ^ _ ^ _ £d + £f j 
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This is not what is usually done in QDVPT, but to pre- 
ceed this way allows to find the direct link to the frozen 
local hole approximation we are looking for. 

In the first step of the FLHA, the model space con- 
figurations and the excited configurations are restricted 
to those determinants which already contain a hole in a 
particular local occupied orbital \<Ph), such that the wave 
operator reduces to 



f?^{h) = \<s> h )(<s> h \+ \Kh) 

a^h,v 



)h\ \ah) 



(16) 



Only single excitations |$^) with respect to the given 
local hole configuration \$h) show up here. Inspection of 
Eq. I|1U|) reveals that these restricted wave operators are 
precisely the leading terms in the full wave operator O, 
in the sense that 



f2 = ^ FLH (^) + n 



where the remainder 

n= E E W» 

a<b,v c^{a,b} 



z\\ab) 



£b 



(17) 



(*d (18) 



rrFLH _ P 
n gh - -*hg 



contains all terms in which the three occupied spin or- 
bitals | a), |6), and |c) are all different, while the leading 
terms are made up by all those contributions in which 
two of the hole indices coincide. 

The final secular matrix of the frozen local hole ap- 
proach, 

h kb := ( ^ LH (4$ ff ] \H\ Cl™(h)[$ h ]) , (19) 
is then given by 

(20) 

+ E W gl w chih (wh\\ch) 

c^h,w 

+ E W g ^ h (ag\\vg) 

+ E E W g ^ chth (ag\\vg)(wh\\ch) 

a^g,v c^h,w 

which is again precisely the leading coinciding-hole- 
indices part of the full effective Hamiltonian (de- 
spite the always present zeroth-order contribution —Fh g ). 
Thus, up to first order perturbation theory the FLHA can 
be understood as a simple neglection of all three-distinct- 
spin-orbital contributions to the wave operator Cl and the 
resulting effective Hamiltonian. 

A further simplification of the FLHA can be achieved 
by treating the (iV-l)-particle system with a frozen lo- 
cal hole on the Hartree-Fock level only. This has already 
been done with some success in our "simplified method" 
to correlation effects in band structures [3(| ■ Orbital re- 
laxation around the localized hole is the only effect which 



can be accounted for in this approach and we want to an- 
alyze here, to which extent the correlation effects in the 
(TV-l)-particle system can be mimicked this way. 

Removal of an electron from a fixed localized occupied 
spin orbital \iph) leads to the following modified Fock 
operator F PLH (/i), 



Fa, 



(ih\\jh) 



(t\F FLH (h)\3) 
with (i\V\j) = (l-SiAFij-iihWjh) 



Sij Si 



(i\V\3) 



(21) 



Up to first order in the perturbation potential V the 
relaxed Hartree-Fock orbitals \(pi(h)) which result from 
diagonalizing F FLH (/i) under the constraint that \(ph) re- 
mains unchanged read 



S 7" Si 



else 



and the corresponding Slater determinant 

(_-i\N-h 
\SCF\ I L ) 



becomes 



(22) 



(23) 



det(^i(/i), . . .,<ph-i(ti),$h+i(ti), ■ ■ .,tpisi(h)) 



i*r> = i**> - £ E \K h 



U\v\ 



0{V 2 



(24) 

Excitations into occupied spin orbitals \ipj) other than 
\tp a ) and \tph) are not possible in the 2hlp configurations 
\^ah)- Thus, j only runs over the virtual orbitals and 
one has 

l^ CF > = |*fc> + E l^)^^ + ^ 2 )(25) 



= fi FLH (/i)|$ h ) + 0(V 2 ) 



(26) 



Obviously, up to first order, the local frozen hole wave 
function obtained on the Hartree-Fock level is identical to 
the fully correlated local hole state \$ h ) = Cl FLH (h)\$ h ) 
of the system (see Eq. (|16[) L This result is somewhat 
surprising, at first glance, because it tells, that what for- 
mally looks like pure orbital relaxation around a frozen 
local defect, is actually bare electron correlation. It also 
tells, that performing a Hartree-Fock calculation around 
a frozen local hole rather than a much more demanding 
wavefunction-based correlation calculations is a rather 
promising approximation. 



V. CONCLUSIONS 

The so-called frozen local hole approximation (FLHA) 
has been analyzed in this study. It presents a method 
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that allows to perform correlation calculations for 
cationic and anionic (iV±l)-electron systems in a much 
more efficient way than running a full quantum chemical 
correlation calculation such as MRCI(SD). In this work 
we were focusing on cationic holes states, but, in a totally 
analogue way, the approximation can also be applied to 
anionic electron attachment states (as done in Ref. |3C|L 
The FLHA is a two step procedure: In the first step, 
the electron to be removed is assumed to reside in a given, 
fixed local occupied orbital and standard wavefunction- 
based correlation calculations are performed in order to 
find the correlation hole around each of these local holes. 
Because of the frozen character of the given hole orbital 
the configuration space is substantially reduced in these 
correlation calculations. Therefore, the approximation 
is a perfect, linear scaling divide-and-conquer approach. 
In the second step, the many-body Hamiltonian matrix 
elements with the correlated local hole states generated 
in the first step are evaluated and the resulting effec- 
tive Hamiltonian matrix is diagonalized. In this way, the 
hole can delocalize over the whole system to form proper 
Bloch states. Hence, the FLHA can be understood as an 
adiabatic treatment, in which the shape of the correla- 



tion hole around each electron hole is kept fixed during 
the final hybridization of the correlated local hole states. 

It has been shown numerically, for two sample systems, 
(H2) n ladder chains and hydrogen terminated linear H- 
(Be)„-H chains, that the FLHA performs astonishingly 
well. Potential energy curves have been calculated on the 
MRCI(SD) level of theory, and deviations of at most 0.1 
eV between a full correlation calculations and the frozen 
local hole approach were found, both, for the van der 
Waals bound H2 chains and for the covalent Be chains. 

Using quasi-degenerate variational perturbation the- 
ory, it was possible to demonstrate that the FLHA in- 
deed assembles all leading terms of a full correlation cal- 
culation. Actually, the approximation consists in a ne- 
glection of all three-distinct-hole-site contributions over 
those where (at least) two hole orbitals coincide. Further- 
more, it could be shown that the correlation calculation 
around a frozen local hole can very well be replaced by 
a numerically much less demanding SCF calculation. In 
first order perturbation theory the resulting orbital re- 
laxations are totally equivalent to the correlation effects 
around a frozen local hole. 
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TABLE I: Analysis of the most relevant configurations in the two lowest correlated cationic hole states l^i) and ^2} of 
Be3fi2. The first columns give the occupation numbers of the involved spatial orbitals (Fig. 0. The indices 5* and T refer to 
spin-adapted configurations with a singlet and triplet-like linear combination of the remaining electrons in a and b, respectively. 
The next two columns show the CI coefficients ai and f3i of the approximate correlated hole states [\t a ) and l^b), respectively, 
as defined in Eq. @, followed by the CI coefficients of the hole states and | \£f ) after the diagonalization step 

0. They are to be compared to the CI coefficients aff ci and afj 01 of the "exact" MRCI wave functions |*^i\ CI } and |*^3, CI >. 
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